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^D ' In this paper we study emergent universe model in the context of a self interacting Jordan-Brans- 

^ . Dicke theory. The model presents a stable past eternal static solution which eventually enters a 

o, 

^' , phase where the stability of this solution is broken leading to an inflationary period. We also 

^D ' establish constraints for the different parameters appearing in our model. 
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I. INTRODUCTION 



Cosmological inflation has become an integral part of the standard model of the universe. 
Apart from being capable of removing the shortcomings of the standard cosmology, it gives 



see 



important clues for large scale structure formation. The scheme of inflation 
for a review) is based on the idea that there was a early phase in which the universe evolved 
through accelerated expansion in a short period of time at high energy scales. During this 
phase, the universe was dominated by a potential t/(^E') of a scalar field \E', which is called 
the inflaton. 

Singularity theorems have been devised that apply in the inflationary context, showing 
that the universe necessarily had a beginning (according to classical and semi-classical the- 
ory) y, [tI, IsI, y, |lO|. In other words, according to these theorems, the quantum gravity era 
cannot be avoided in the past even if inflation takesplace. However, recently, models that 
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escape this conclusion has been studied in Refs. [ll|, 
models do not satisfy the geometrical assumptions of these theorems. Specifically, the the- 
orems assume that either i) the universe has open space sections, implying A; = or — 1, 
or ii) the Hubble expansion rate H = a/ a is bounded away from zero in the pas t, H > 0, 



where a is the scale factor. In particular, Refs 
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closed models in which k = +1 and H can become zero, so that both assumptions i) and ii 
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of the inflationary singularity theorems are violated. The models studied in Refs. 
obey general relativity, contain only ordinary matter, and (minimally coupled) scalar fields. 
In these models, the universe is positively curved and initially in a past eternal classical 
Einstein static (ES) state that eventually evolves into a subsequent inflationary phase. Such 
models are appealing since they provide specific examples of nonsingular (geodesically com- 
plete) infiationary universes. Furthermore, it has been proposed that entropy considerations 
favor the ES state as the initial state for our universe |l9l . |20(] . 

However, the models based on general relativity with ordinary matter suffer from a num- 
ber of important shortcomings. In particular, the instability of the ES state (represented 



by a saddle equilibrium point in the phase space of the system, see [13|, Il6|, Il7|, Il8| ) makes 
it extremely difficult to maintain such a state for an infinitely long time in the presence of 
fiuctuations, such as quantum fiuctuations, that will inevitably arise. As in the emergent 
universe scenario, it is assumed that the initial conditions are specified such that the static 



configuration represents the past eternal state of tlie universe, out of which the universe 
slowly evolves into an inflationary phase. The instability of the ES solution ensures that 
any perturbations, no matter how small, rapidly force the universe away from the static 
state, thereby aborting the scenario. Some models have been proposed to solve the sta- 
bility problem of the asymptotic static solution. They consider non-perturbative quantum 
corrections of the Einstein field equations, either coming from a 'semiclassical' state in the 
framework of loop quantum gravity (LQG) |13l . Il7l | or braneworld cosmology with a timelike 



extra dimension 
matter pj,ll5j. 



18| . Other possibilities to consider are the Starobinsky model or exotic 



The Jordan-Brans-Dicke (JBD) 2l|] theory is a class of models in which the effective 
gravitational coupling evolves with time. The strength of this coupling is determined by 
a scalar field, the so-called Brans-Dicke field, which tends to the value G~^, the inverse of 
the Newton's constant. The origin of Brans-Dicke theory is in Mach's principle according 
to which the property of inertia of material bodies arises from their interactions with the 
matter distributed in the universe. In modern context, Brans-Dicke theory appears naturally 
in super gravity models, Kaluza-Klein theories and in all the known effective string actions 
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In this work we consider a JBD model and determine whether such a model could fit 
the general characteristics of an emergent universe scenario: a stable static past asymptotic 
solution followed by a period of de Sitter inflation. Here, instead of employing a time-like 
extra dimension or examining a past eternal static solution that lies in the semiclassical 
quantum gravity regimen of the theory, we work more conventionally, keeping our model 
just at the classical level, in the spirit of Refs. 
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The paper is organized as follows. In Sect. [TTlwe review briefly the cosmological equations 
of the JBD model. In Sect. Illll the existence and nature of a static solution is discussed. In 
Sect. IIVI we study the dynamics that lead to the emergence of an inflationary universe. In 
Sect. |V]we present a specific model that satisfies the requirements of an emergent universe 
in the scheme of JBD theories. In Sect. I VI I we summarize our results. 



II. THE MODEL 



We consider the following JBD action 2l|] for a self-interacting potential and matter, 
given by 



S = d^x \/—g 



1 ^ ^ 1 w 
- ^R 

2 2 $ 



$ i? - - - 9^$ a^$ - K($) + c{m) 



(1) 



where >C(^I/) denote the Lagrangian density of the matter 

i? is the Ricci scalar curvature, $ is the JBD scalar field, w is the JBD parameter and 
y($) = V is the potential associated to the field $. Here \1/ is the standard infiaton field 
and U{^) its effective potential. In this theory 1/$ plays the role of the gravitational 
constant, which changes with time. This action also matches the low energy string action 

for w = -1 y. 

The Friedmann-Robertson- Walker metric is described by 

ds^ = df - a{tf dQl , (2) 

where a{t) is the scale factor, t represents the cosmic time and dQ^. is the spatial line element 
corresponding to the hypersurfaces of homogeneity, which could represent a three-sphere, a 
three-plane or a three-hyperboloid, with values k=l, 0,-1, respectively. From now on, we 
will restrict ourselves to the case k = 1. 

Using the metric ([2]), with k=l, in the action ([1]), we obtain the following field equations: 

a . 1 $ ^6 w f^Y V P 

a a? (^ $2\$/$<l> ^ ' 

$ + 3/7$ = [P ~ ^^l + -J— \2V-^V'] , (5) 

and the conservation of energy-momentum implies that 

p + 3H{p + P) = 0, (6) 



or equivalently 

i + 3M = -^, (7) 

where V = dV{^)/d^. Dots mean derivatives with respect to time, units are such that 
Svr G = 1 and c = h = 1. 

Here the energy density p and pressure p are given by 

p=^ + f/(vl/), (8) 

and 

P = — - t/(vl/) . (9) 

We could write an effective state equation for this scalar field given by P = (7 — 1) p, 
where the equation of state "parameter", 7, is defined by 
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In the situation where the scalar potential of the inflaton field is constant, the state 
parameter is a function only of the scale factor. 

III. STATIC UNIVERSE 

A. Static solution in JBD model 



In the JBD model the static universe is characterized by the conditions a = a^ = Const., 
do = = ao and $ = $0 = Cte., $0 = = $o- Following the same scheme as the 
static Einstein model, we are going to consider that the matter potential f/(\E') becomes 
asymptotically flat in the limit \E' -^ — cxd, that is f/(^E') = Uq = Const.. In this limit, the 
initial conditions are specified such that the static configuration represents the past eternal 
state of the universe, out of which the universe slowly evolves into an inflationary phase. 
Then from equations ([3]) to ([5]) and using the equation of state P = (7 — 1) p, we obtain the 
following equations of stability 



= ^!L + J:1_1 

3$o 3<l>o al' 



= (4-37o)po + 4Vo-2<l'oK)', (H) 



^ , f IN Po V'c 







= - + (70 - 1) ^ ^ 



where Vo = V{%), V^ = {dV{^)/d^)^=^, and 70 = 2 ('l - 

The equations (fTTj) are satisfied if the following conditions are fulfilled 
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(12) 
(13) 

Po = V^<l>o-Vo. (14) 

We can obtain the velocity at which the scalar field \E' is rolling along the constant 
potential Uq as a function of the static values of the scale factor ao, Brans-Dicke field $0 
and energy density po 

^ = 2^. (15) 

Note that in order to obtain a static solution we need to have a non-zero JBD potential 
with a non- vanishing derivative at the static point $ = $o- The original Brans-Dicke model 
corresponds to V{^) = 0. However, non-zero V^($) is better motivated and appears in many 
particle physics models. In particular, V{^) can be chosen in such a way that $ is forced 
to settle down to a non-zero expectation value, $ -^ mi, where nip = lO^^GeV is the value 
of the Planck mass today. On the other hand, if V{^) fixes the field $ to a non- zero value, 
then time-delay experiments place no constraints on the Brans-Dicke parameter w [29(|. 



B. Oscillations 

As we have mentioned above, one important point that we have to determine is whether 
the static JBD solution found in the previous section corresponds to an stable solution. In 
order to see this, let us consider small perturbation about the static solution for the scale 
factor and the JBD field. We set 



a(t) = ao [1 + e{t)] 



(16) 



and 



m = %[i+pm, (17) 

where e -C 1 and /3 -C 1 are small perturbations. By introducing the expressions (IT6|l and 
( TT7|) into Eq. (jl]) and Eq. ([5]) and retaining only at the linear order in e and (3 we obtain the 
following coupled equations: 
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where Vq" = {d^V{<l>)/d<l>^)^=^,. 










Here, we have used that 











(18) 



(19) 



p = po + 5p{e) ^ po - 37o Po £ ■ 



(20) 



and 



7 = 7o + hi^) ~ 7o - 6 e . 

Po 



(21) 



From the system of Eqs. OlSI) and (TT9l) we can obtain the frequencies of small oscillation 



u. 



al{3 + 2w) 



al%V^' - 2(2 + 3w) 



± ^J[al%V^']^ + 4:al%V^'{3 + 2w) + 8wi3 + 2w) . (22) 

The static solution is stable if tu^ > 0. Assuming that (3 + 2w) > 0, we found that the 
following inequalities must be satisfied in order to have a stable static solution 



o<a2<i.ov;;'<|, 



and 



< w < 



(23) 
(24) 



/9 - 6a2$oVo" + (3 + a^$oK") 

These inequalities restrict the parameters of the model. The first imposes a condition on 
the JBD potential, specifically for its first and second derivatives: < Vq < Vq/{2^q). The 
second inequality restrictions the values of the JBD parameter. Notice that this inequality 
imposes that w < 0. JDB models with negative values of^u' have been considered in the 
context of late acceleration expansion of the universe 
limits of string theory 28 1. 
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3l| , but also appear in low energy 



In our case we are going to choose the JBD potential in such a way that $ will be forced 
to stabilize at a constant value $/ at the end of the inflationary period, see next section. 
Then, we can recover General Relativity by setting $/• = m'i, therefore whatever we choose 



for w in our model, it does not contradict the solar system bound on w 291.132 1. 



IV. LEAVING THE STATIC REGIME 

The discussion of the previous section determined the behavior of the model in the case of 
a constant potential for the scalar field \1/. However, any realistic inflationary model clearly 
requires the potential to vary as the scalar field evolves. Here, following Ref. 13| and with 
the emergent inflationary models in mind, we consider a general class of potentials that 
approach a constant ?7o as \I^ ^ — oo and rise monotonically once the value of the scalar 
field exceeds a certain value. 

The overall effect of increasing the potential is to distort the equilibrium behavior ex- 
pressed by Eqs. f lT2|) and flTSj) . The inclusion of the derivative term in the equation of the 
scalar field ^ produce changes in its equilibrium velocity, Eq. (TTSll . breaking the static solu- 
tion. In particular, the field \1/ decelerates as it moves further up the potential, subsequently 
reaching a point of maximal displacement and then rolling back down. If the potential has 
a suitable form in this region, slow-roll inflation will occur. 

On the other hand, in the slow-roll regimen the scalar potential evolves slowly. In that 
case we can consider f/(\E') ~ const. = Umf. Then, as mentioned in Ref. [33|, Eqs. ([3]) and 
([5]) have an exact static solution for a particular value of $, driving a de Sitter expansion. 
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This occurs when the right hand side of Eq. ([5]) becomes zero. Denominating this quasi-static 
value of the JBD field as $, it satisfies the following condition 

AU,nf + ^V{^)-2^V'{^)=0. (25) 

Then, once the scalar field starts to move in the slow-roll regime, the JBD field goes to 
the value $ and the universe begins a de Sitter expansion with 



3$ 



U,nf + ^($)J . (26) 

For example, the JBD potential V{^) = A(<l'^ — z/)^ satisfies this condition, see Ref. 33|. 
In the next section we consider a specific potential which satisfies Eq. (1251) together with the 
requirements of a static solution described in Sect. IIIII 

Finally, during the evolution of \1' over f/(\E') to zero, the JBD field evolves slowly to its 
final value $/, at which the expression 2V — (pV vanishes. We consider the value $/ as the 
current value of the JBD field. 

Also, we can determinate the existence of the inflationary period by introduce the dimen- 
sionless slow-roll parameter 



e = — - 



V + U 



H^ (3 + 2w) ' 2\V + Uj ' {3 + 2w){V + U) [V + U "J" ^'^^^ 

Then, the inflationary regime takes place if the parameter e satisfies the inequality e < 1, 
a condition analogous to the requirement a > 0. We note that if Eq. (I25l) is satisfied (i.e. 
$ = $) and the scalar potential f/(^l') satisfies the requirement of an inflationary potential, 
we get e < 1. In the next Section we study a particular model that follows the behavior 
described above. 

V. A SPECIFIC MODEL OF AN EMERGING UNIVERSE 

From a dynamical point of view, the emergent universe scenario can be realized in the 
context of JBD if the scalar potential t/(\E') satisfies a number of weak constraints. Asymp- 
totically, it should have a horizontal branch as \E' ^ — oo such that dU /d^ — > and increase 
monotonically in the region ^ > '^ grow where without loss of generality we may choose 
^ grow = 0. The reheating of the universe imposes a further constraint. There should be a 



global minimum in the potential at Umin = 0, if reheating is to proceed through coherent 
oscillations of the infiaton. The region of the potential that drove the inflationary expansion 
is then constrained by cosmological observations, as in the standard scenario. 

Motivated by the former discussion, we consider the following potential as an example: 



U{^) = Uo [exp i^P^/Vsj - ij , 
which exhibits the generic properties described above, see Fig. ([1]). 



(28) 



x10 



12 



U(^) 




FIG. 1: The plot shows an emergent scalar potential that allows for conventional reheating. 



Potentials of this form have been considered previously, not only in the context of the 



emergent universe 



11 



fl^ 



Il3l | but in a number of different settings, including cases that 
introduce higher-order curvature invariants into the Einstein-Hilbert action. Such correc- 
tions are required when attempting to renormalize theories of quantum gravity 34|- They 
also arise in low-energy limits of superstring theories 35|. In general, these theories are 
conformally equivalent to Einstein gravity plus a minimally coupled, self-interacting scalar 
field. In particular, potentials with the structure of Eq. ( I28l) can be obtained from theories 

12]. In general, all these 



that include a R^ term in the action, where R is the Ricci scalar 



potentials possess a global minimum at ^ = 0. Following Ref. [l3|], we take Ui 
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and /3 = 0.1 as representing typical values satisfying the constraints imposed by the WMAP 
satellite [Sy, |37| • 

As an example of a Brans-Dicke potential that satisfies the condition of static solution, 
Eqs. f l25]) and i^^, we consider the following polynomial potential 

V{<^) = K) + ^ ($ - $o) + ^ 5 ($ - $o)' + ^ C ($ - $o)' , (29) 

where $o correspond to the value of the JBD field at the static solution. 

Following the discussion of Sect. llVt we choose the parameter C in order to force $ to 
settle down to a non-zero expectation value, $ — > $j. Then we have: 



C = -12 



2Vo - 2A(!>o + B ^l + A^f - B % "^f 



The parameters Vo,A,B are fixed in order to obtain a static solution at ^ = ^o, a = ao 
and p = Po 



Vo = Po , (30) 

ao 

>4 = |, (31) 

where we have introduced the dimensionless parameter x, satisfying < x < 3/2. The JBD 
parameter requires: 



^<w <-^./3^2^-\i3 + x). (33) 



The static energy density is given by: 



Po = ^ + Uo. (34) 



% 



In order to obtain a numerical solution we take the following values for the parameters 
in the JBD potential $o = 0.9, $/ = 1, ao = 5.4 x 10^, x = 1 ^ind w = —1.45, where we 
have used units in which SttG = 1. These particular parameters satisfy all the constraints 
discussed previously. On the other hand, in order to consider the model just at the classical 
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FIG. 2: Behavior of the scale factor a, the JBD field $ and the scalar field ^ as function of the 
cosmological time t. 

level we have to be out of the Planck era. This imposes the following conditions: p < $(t)^ 
and V{^) < $(t)^, which are satisfied with the values of the parameters mentioned above. 
Now, let us consider a numerical solution corresponding to a universe starting from an 
initial state close to the static solution. The whole evolution of the model is shown in 
Fig. ([2]), where we can notice that during the part of the process where the scalar field \E' 
moves in the asymptotically flat region of the potential f/(\I/), the universe remains static, 
due to the form of the JBD potential, Eq. (129|) . This means that the scale factor and the 
JBD field do not evolve with the cosmological time keeping their equilibrium values Oq and 
$0 during this period. 
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FIG. 3: Plot showing a, ^ and ^ near the beginning of the inflationary period. 

The static regimen finishes when the scalar field moves past the minimum of its potential 
(near the value \E' ~ 0) and begins to decelerate as it moves up its potential. During this 
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period the static equilibrium is broken, and the scale factor and with the JBD field start 
to evolve. Finally, the scalar field starts to go down the potential f/(\l') in the slow-roll 
regimen. The details of the last part of this process is shown in Fig ([3]), where we note that 
at the moment when the scalar field starts to roll back down its potential, the JBD field 
attains its quasi-static value $ ~ 1, and the scale factor starts its infiationary expansion. 

On the other hand, a numerical solution corresponding to a universe starting from an 
initial state not in the static solution but close to it, presents small oscillations around the 
equilibrium values, as shown in Fig (j4]). This tells us that the static solution is stable. 
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FIG. 4: Graph showing the scale factor and the JBD field oscillating around their equilibrium 
values ao and <1>0) while the scalar field follows the constant potential. 



VI. CONCLUSIONS 



In this paper we have studied a Jordan-Brans-Dicke model and we have determined 
whether that model could display the general characteristics required for an emergent uni- 
verse scenario. That is a stable static past asymptotic solution followed by a period of de 
Sitter infiation. 

The original idea of an emergent universe |ll|] is a simple closed infiationary model in 
which the universe emerges from an Einstein static state with radius Oq >> Lp, infiates and 
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is then subsumed into a hot Big Bang era. The attractiveness of the proposed model is that 
one can avoid an initial quantum-gravity stage if the static radius is larger than the Planck 
length. However, this model suffers from the problem of instability of the Einstein static 



n 



state (see Refs. [l3|, ll6|, [l7|, [l8|) which makes it extremely difficult to maintain its state for 



an infinitely long time in the presence of fiuctuations, such as quantum fiuctuations, thereby 
aborting the scenario. 

In this work, we have provided an explicit construction of an emergent universe scenario, 
which presents a stable past eternal static solution and brings us the possibility of avoiding 
an initial quantum-gravity stage if we chose the static radius to be larger than the Planck 
length. 

In particular, we have considered a JBD theory with a self interacting potential and 
matter content corresponding to a scalar field. 

In the first part of the paper, we studied static solutions. In order to do so we determined 
the characteristics of the JBD potential and we took a constant scalar matter potential. In 
particular, we have found a static solution in which the JBD potential had a non-zero value 
and a positive derivative at the static point $ = $o- In determining the stability of this 
solution, we have calculated the real frequencies of small oscillation about the static solution. 
This imposed a bound for the first and second derivative of the JBD potential at the static 
point < V^" < V"q'/(2$o)- A restriction on the value of the JBD parameter w was also 
obtained (see Eq. ( |24l) ). In this way, we have shown that it is possible to obtain a past 
eternal universe in a JBD model, depending upon the characteristics of the JBD potential 
and the Brans-Dicke parameter w. 

In the second part of the paper, we studied the possibility that our model present a 
past eternal static solution, which eventually enters a phase where the stability of the static 
solution is broken by changing the matter scalar field potential, thereby leading to a phase 
of infiation. 

In our model, the mechanism that enables the universe to emerge depends on the form 
of both potentials: the JBD potential V^(0) and the matter scalar field potential, U{^). For 
the scalar field potential it is required that it asymptotically approaches a constant value as 
\1/ -^ — oo, and in order to break the cycles it should grow in magnitude for larger \1/. The 
JBD potential must satisfy similar requirements to that described in Ref. |33|. 

In the third part of the paper, we studied a particular matter scalar potential, similar 
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to the one used in the context of emergent universe [ll|, Il2|, Il3|, and a polynomial JBD 
potential, which satisfies the requirement of static stable past eternal solution followed by a 
period of infiation, Eq (1^ . 

We obtained numerical solutions for a universe starting from an initial state close to 
the static solution. The numerical solutions showed a behavior just like that discussed in 
previous sections. In particular we have found that when the scalar field \1/ moves in the 
asymptotically fiat region of the potential ^(^E^), the scale factor and the JBD field experience 
small oscillations about their equilibrium points. After that, when the scalar field passes 
the minimum of its potential and begins to decelerate as it moves up its potential, we found 
that the static equilibrium is broken, and the scale factor and JBD field start to evolve. 
In particular, the numerical solution shows that when the scalar field starts to go down 
the potential f/(^I/), the JBD field gets its quasi-static value and the scale factor begins a 
quasi-exponential expansion. 

We should note that a more detailed analysis of this process could be done by using a 
dynamical system approach. We expect to return to this point in the near future. 
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